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Abstract. 

We use OpenGGCM global MHD simulations to study the nightside magnetospheric, 
magnetotail, and ionospheric responses to interplanetary (IP) fast forward shocks. Three 
cases are presented in this study: two inclined oblique shocks, hereafter IOS-1 and lOS- 
2, where the latter has a Mach number twice stronger than the former. Both shocks have 
impact angles of 30° in relation to the Sun-Earth line. Lastly, we choose a frontal per¬ 
pendicular shock, EPS, whose shock normal is along the Sun-Earth line, with the same 
Mach number as IOS-1. We find that, in the IOS-1 case, due to the north-south asym¬ 
metry, the magnetotail is deflected southward, leading to a mild compression. The ge¬ 
omagnetic activity observed in the nightside ionosphere is then weak. On the other hand, 
in the head-on case, the EPS compresses the magnetotail from both sides symmetrically. 

This compression triggers a substorm allowing a larger amount of stored energy in the 
magnetotail to be released to the nightside ionosphere, resulting in stronger geomagnetic 
activity. By comparing IOS-2 and EPS, we find that, despite the IOS-2 having a larger 
Mach number, the EPS leads to a larger geomagnetic response in the nightside ionosphere. 

As a result, we conclude that IP shocks with similar upstream conditions, such as mag¬ 
netic field, speed, density, and Mach number, can have different geo effectiveness, depend¬ 
ing on their shock normal orientation. 


1. Introduction 

Interplanetary (IP) shocks are ubiquitous features of the 
solar wind [Burlaga, 1971]. As they encounter Earth they 
interact with the magnetosphere, causing disturbances that 
can be seen everywhere in the magnetosphere. Because the 
disturbances alter the magnetospheric current systems, the 
ionosphere is also affected, and the magnetic field on the 
ground is perturbed as well. The most dramatic shock- 
induced ground perturbations are the storm sudden com¬ 
mencements (SSCs), which are driven by a strong IP shock 
preceding a geomagnetic storm driven by coronal mass ejec¬ 
tions (CMEs). 

The interaction of IP shocks with the Earth’s magneto¬ 
sphere is both complex and important. For example, the 
shock-shock interactions such as between an IP shock and 
the Earth’s bow shock may occur in many contexts, for 
example, in the heliosphere and in astrophysical systems. 
Such remote interactions are difficult to observe, but can be 
readily observed with in-situ measurements in the magne¬ 
tosphere. On the other hand, strong IP shock impacts on 
the magnetosphere have substantial space weather effects, 
for example, they produce geomagnetically induced currents 
(GICs), which can impact power grids [Bolduc, 2002; Kap- 
penman, 2010], and they can energize particles in the in¬ 
ner magnetosphere [Hudson et al., 1997; Zong et al, 2009]. 
Echer et al. [2004] reported that 22% of all interplanetary 
shocks are intensely geoeffective while 35% are moderately 
geoeffective. Thus, the study of IP shock impacts is of fun¬ 
damental and also of practical importance. 

At 1 AU (astronomical unit) IP shocks are almost exclu¬ 
sively fast shocks. Slow shocks may exist at closer distances 
to the Sun, but they are subject to Landau damping [Richter 
et al., 1985]. IP shocks may further be classified as propa¬ 
gating away from the Sun, i.e., forward shocks, or as propa¬ 
gating towards the Sun, i.e., so-called reverse shocks. Since 
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the solar wind speed is almost always supermagnetosonic, 
a reverse shock will still propagate away from the Sun in 
the Earth’s frame. IP shocks may then be further classi¬ 
fied by their strength, i.e., their Mach number in the solar 
wind reference frame, their orientation, and by the orienta¬ 
tion of their upstream field with respect to the shock nor¬ 
mal. The compression ratio, i.e., the ratio of downstream 
to upstream plasma density, is an alternative measure for 
the shock strength, which is often more convenient to use 
than the Mach number. Assuming a 7 (specific heat ratio) 
of 5/3 appropriate for a monoatomic gas, the compression 
ratio must lie between 1 and 4 [Priest, 1981]. IP shocks are 
typically weak (compared to planetary bow shocks) with a 
compression ratio between 1.2 and 2 [Berdichevsky et al., 
2000]. The Rankine-Hugoniot relations provide the jump 
conditions between the upstream and downstream plasma 
parameters in the MHD context [Jeffrey and Taniuti, 1964; 
Priest, 1981; Parks, 2004]. Thus, given the upstream plasma 
and field parameters, as well as the shock normal, the down¬ 
stream parameters can easily be calculated. The inverse cal¬ 
culation, i.e., determining the shock speed and orientation 
from measured upstream and downstream values, is gen¬ 
erally a much more difficult problem, because the critical 
parameters that cannot be measured directly, such as the 
shock normal, depend in a very sensitive manner on the 
upstream and downstream plasma and field measurements. 
Near 1 AU, IP shocks can generally be assumed to be planar 
structures on the scale size of the Earth’s magnetosphere. 
For example, Russell et al. [1983] found the assumption of 
planarity consistent with measurements from four widely 
spaced solar wind monitors. Also using the assumption of 
shock planarity, Russell et al. [2000] estimated the shock 
normal orientation of a large IP shock with accuracy by 
comparing results of different IP shock normal determina¬ 
tion methods. 

The shock normal vector determines how the shock prop¬ 
agates through the heliosphere. Normals of most IP shocks 
generated by CMEs at 1 AU are concentrated near the Sun- 
Earth line [Richter et al, 1985]. However, shocks driven 
by corotating interaction regions (CIRs), as a result of the 
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slow solar wind compression by a fast stream, have normals 
inclined in relation to the Sun-Earth line [see, e.g., Siscoe, 
1976; Pizzo, 1991, and references therein]. For CIR-driven 
shocks, the normal angles in the azimuthal direction in re¬ 
lation to the solar coordinate system are generally equal or 
larger than the inclination angle [Siscoe, 1976; Pizzo, 1991]. 

Here, we are primarily concerned with the effect that the 
orientation of an IP shock has on its geoeffectiveness. Sim¬ 
ilar studies have been done in the past, but none have ad¬ 
dressed the parameter space that we are concerned with. 
For example, Jurac et al. [2002] analyzed more than 100 
IP shocks in terms of the angle between the shock nor¬ 
mal and the upstream interplanetary magnetic field (IMF), 
also called obliquity. They found that there is about a 40% 
chance for an IP shock to be followed by an intense mag¬ 
netic storm when the obliquity is 90°. These shocks are more 
likely to result in an intense storm than shocks whose nor¬ 
mal are oblique to the IMF direction. Takeuchi et al. [2002] 
showed that the shock orientation plays an important role 
in controlling the rise time duration of an SSC, such that 
the rise time is longer in the case of larger impact angles. 
The unusual SSC rise time observed, ~30 minutes, was ex¬ 
plained by a gradual compression observed in the magne¬ 
tosphere due to the fact that the shock normal had to be 
largely inclined to the dusk-ward direction. The interaction 
of IP shocks inclined in relation to the Sun-Earth line with 
the bow shock was also investigated by Grib and Pushkar 
[2006]. They solved the Rankine-Hugoniot conditions nu¬ 
merically for different shock inclinations and one of their 
most important results was that, for shock normal inclina¬ 
tions between 60° and -60°, the density changed from dusk 
to dawn in the bow shock in the case where the discontinu¬ 
ity was a fast forward shock (FFS). Using ACE and Wind 
satellite data from 1995 to 2004, Wang et al. [2006] reported 
that, in a survey of nearly 300 FFSs, 75% of them were fol¬ 
lowed by SSCs observed on the ground. They also found that 
the shock impact angle plays an important role in determin¬ 
ing the SSC rise time, as previously suggested by Takeuchi 
et al. [2002]. When the shock speed (shock strength) was 
hxed, the more parallel the shock normal with the x-line, the 
smaller the rise time. The same occurred when they fixed 
the shock inclination and changed the shock speed. The 
faster the shock, the shorter the SSC rise time. 

Several authors have studied the interaction of IP shocks 
with the Earth’s magnetosphere in the context of numerical 
MHD simulations. However, almost always the IP shock hit 
the bow shock at the subsolar point head-on. For instance, 
Ridley et al. [2006] simulated an extreme IP shock driven 
by a Carrington-like CME [Manchester et al., 2006] that 
pushed the magnetopause toward the Earth to the limit of 
their code boundary, which was at 2 Re . They also observed 
a secondary shock wave reflected back by the magnetopause 
that encountered the bow shock that was moving inward. 
Then the combined motion propagated down the flanks of 
the magnetosphere. Similar bow shock Earth-ward motion 
was observed by Safrdnkovd et al. [2007] in their numerical 
simulation for a much weaker IP shock as well. The interac¬ 
tion of IP shocks with the Earth’s magnetosphere can also 
lead to generation of two ionospheric current systems [Ridley 
et al., 2006; Guo and Hu, 2007; Samsonov et al, 2010]. The 
appearance of an anomalous region I current, which flowed 
oppositely to the region I current, followed by a frontal IP 
shock impact with no IMF was found by Guo and Hu 
[2007]. This anomalous region I current formed at noon, 
developed and then moved toward the evening side until it 
vanished. Such a response depends on the strength of the 
IP shock. In other MHD simulations concerning the mag¬ 
netosheath, three new discontinuities appeared downstream 
from the bow shock in addition to the FFS: a forward slow 
expansion wave, a contact discontinuity, and a reverse slow 
shock [Koval et al., 2006; Samsonov et al., 2006, 2007]. Sam¬ 
sonov et al. [2010] simulated the interaction of an IP shock 


with the magnetosphere in an artihcial case with a north¬ 
ward IMF. The IP shock normal was aligned with the Sun- 
Earth line. They observed an intensification of two iono¬ 
spheric current systems (similar to the preliminary and main 
impulse currents) that coincided in time with the intensih- 
cation of two corresponding magnetospheric dynamos. 

Guo et al. [2005] performed a global numerical MHD sim¬ 
ulation with different shock normal orientations to study the 
interaction of IP shocks with the Earth’s magnetosphere. 
They simulated two different cases with a Parker-spiral IMF 
orientation with no Bz component. In their first case, the 
shock normal was parallel to the Sun-Earth line, and in their 
second case the shock normal was oblique to this line with 
an angle of 60°. They found that the inclined shock nor¬ 
mal led to a longer evolution time of the system. Although 
both systems evolved from the same initial conditions in 
their numerical simulations, as can be seen in their Figure 
4, they did not find any signihcant difference in the hnal 
quasi-steady state of the systems. Similar results were also 
found by Wang et al. [2005]. More recently, similar results 
have been found by Samsonov [2011] as well. He presented a 
solution and analysis of the problem of an inclined IP shock 
incident on and propagating through the Earth’s magne¬ 
tosheath. He showed that inclined IP shocks with normals 
in the equatorial plane result in a dawn-dusk asymmetry in 
the magnetosheath and predicted that this effect should be 
present in observations of sudden impulse inside the magne¬ 
tosphere and on the ground. 

Here, we show that the IP shock normal orientation is 
a critical parameter determining the geoeffectiveness of IP 
shocks. Preliminary results of this work were presented by 
Oliveira and Raeder [2013]. IP shocks that have their nor¬ 
mal aligned with the Sun-Earth line have the largest geoef¬ 
fectiveness, because they compress the magnetosphere from 
all sides at the same time. By contrast, if the shock nor¬ 
mal makes a large angle with the Sun-Earth line in the x-z 
plane in GSE (Geocentric Solar Ecliptic) coordinate system, 
the north-south asymmetry of the impact pushes the plasma 
sheet either to the north or to the south without much com¬ 
pression, which leads to a much weaker response. 

In the following, we describe the model in section 2, in 
section 3 we present our results, and in section 4 we sum¬ 
marize and discuss our results. 

2. OpenGGCM Model 

We use the Open Geospace General Circulation Model 
(OpenGGCM) to study the impact of IP shocks on the mag¬ 
netosphere. The OpenGGCM is a global coupled model of 
the Earth’s magnetosphere, ionosphere, and thermosphere. 
It is available at the Community Coordinated Modeling Cen¬ 
ter (ccmc.nasa.gov) as a community model for model runs on 
demand [Rastdtter et al, 2013; Pulkkinen et al, 2011, 2013]. 
The magnetosphere part solves the MHD equations as an 
initial-boundary-value problem. The MHD equations are 
solved to within -^3 Re of Earth. The region within 3 
Re is treated as a magnetosphere-ionosphere (MI) coupling 
region [Oliveira, 2014] where physical processes that cou¬ 
ple the magnetosphere to the ionosphere-thermosphere sys¬ 
tem are parameterized using simple models and relation¬ 
ships. The ionosphere-thermosphere system is modeled us¬ 
ing the NOAA CTIM (Coupled Thermosphere Ionosphere 
Model [Fuller-Rowell et al., 1996; Raeder et al., 2001a]). The 
OpenGGCM has been described with some detail in the lit¬ 
erature [see, e.g. Raeder et al, 2001b; Raeder, 2003; Raeder 
et al, 2008]; we thus refer the reader to these papers. The 
OpenGGCM has been used for numerous studies, including 
studies of substorms [Raeder et al., 2001b; Ge et al., 2011; 
Gilson et al., 2012; Raeder et al., 2013], storms [Raeder et al., 
2001a; Raeder and Lu, 2005; Rastdtter et al., 2013; Pulkki¬ 
nen et al., 2013], and, most relevant for this study, for the 
study of IP shock impacts [Shi et al, 2013]. 
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3. Shock Impacts 

3.1. Simulation setup 

We use exclusively GSE coordinates in our simulation in¬ 
put data. The numerical box extends from the Earth 30 
Re in the Sun-ward direction, and 300 Re down the tail. 
In the directions perpendicular to the Sun-Earth line, i.e., 
in the y and z directions, the numerical box extends to ±50 
Re. The numerical grid is non-equidistant Cartesian and is 
divided into 610x256x256 grid cells, such that the highest 
resolution is closest to Earth. Specifically, the grid resolu¬ 
tion is 0.15 Re within a distance of 10 Re radially from 
the Earth. The inner boundary, where the magnetosphere 
variables connect via field-line mapping to the ionosphere is 
located at 3 Re. 

In this paper we only consider IP shocks for which the 
shock normal lies in the GSE x-z plane. Furthermore, we 
assume that the IMF also has no GSE y-component. Thus, 
the shock geometry relative to the Earth and to the IMF 
depends exclusively on two angles. First, depending on the 
shock normal relative to the upstream (relative to the shock) 
magnetic field direction, a shock can be classified as per¬ 
pendicular, oblique, or parallel [Burlaga, 1971; Tsurutani 
et al, 2011]. As is often found in the literature [Burlaga, 
1971; Tsurutani et al., 2011], when 0° < 9b„ < 30°, the 
shock is classified as almost parallel. In the cases in which 
30° < 9b„ < 60°, the shock is said to be oblique. Finally, 
when 60° < 9b„ < 90°, the shock is classified as almost per¬ 
pendicular. In particular for this paper, the shock is named 
perpendicular when 9b„ ~ 90°. Second, in relation to the 
Earth’s system of reference, the shock normal is decomposed 
in terms of two angles: the angle 9x„ between the shock nor¬ 
mal and the Sun-Earth line, and the angle in the x-y 
plane that completes the set, following the notation of Vinas 
and Scudder ]1986]. For our simulations presented here, ifiy„ 
is always 90°. We also assume average solar wind conditions, 
with a particle number density of 5 cm“®, thermal plasma 
pressure of 20 pPa, magnetic field magnitude of 7 nT, and 
background speed of v\ = (-400,0,0) km/s in all cases up¬ 
stream of the IP shocks. We specify the shock strength by 
its compression ratio. As reported by Echer et al. (2003], 
most shocks near Earth have a compression ratio of the or¬ 
der of 2 during solar maximum. Here, we then choose a 
compression ratio value of 1.5 in order to have a mild shock. 
The MHD code input is set as follows. We transform the up¬ 
stream initial conditions into the shock frame, calculate the 
downstream parameters, and subsequently transform them 
back to the Earth’s system of reference. We ran simula¬ 
tions of different FFSs with different shock normal orienta¬ 
tions, obliquities, shock speeds, Mach numbers, and IMF Bz 
pointing either northward or southward. For brevity, we se¬ 
lect three FFSs with different shock normal orientations and 
Mach numbers to discuss in detail. In the first case, the IP 
shock has its normal inclined with an angle 9x„ of 30° with 
the GSE x-axis toward the south, 9b„ =51°, shock speed 
Va — 380 km/s, and Mach number of 3.7, i.e., an inclined 
oblique shock, hereafter IOS-1. The second case corresponds 
to an FES, also inclined and oblique here called IOS-2, with 
9x„ = 30°, 9b„ = 45°, Va = 650 km/s, and Mach number 
of 7.4. In the last case the shock normal has 9x„ = 0° and 
was perpendicular to the IMF, i.e., a frontal perpendicu¬ 
lar shock (FPS), with Va = 650 km/s, and Mach number 
of 3.7. These shock speeds are consistent with the obser¬ 
vations reported by Berdiehevsky et al. (2000], where most 
FFSs have speeds in the range 50-200 km/s in the shock 
frame of reference. Tables 1,2, and 3 show the upstream (1) 
and downstream (2) and other important parameters for the 
three FFSs. The first and second shocks impact the magne¬ 
tosphere (first contact with the magnetopause) at t = 16.45 


minutes, and the FPS reaches the subsolar magnetopause 
at t=18.28 minutes. We also simulated shocks with north¬ 
ward IMF and otherwise identical solar wind conditions. We 
found that the results were similar to the southward cases, 
but with weaker magnetosphere response, with the excep¬ 
tion that transient northward (NBZ) currents occurred 
within a few seconds after the shock impact when it was 
frontal. This effect was already reported by Samsonov et al. 
(2010]. We will thus focus on the case of southward IMF. 

3.2. Results 

As the IP shock impacts the magnetopause, it launches 
waves into the magnetosphere. The phase speed of these 
waves (both Alfven and magnetosonic waves) is generally 
much larger in the magnetosphere than in the solar wind 
and the magnetosheath. Therefore, the amplitude of the 
waves diminishes in the magnetosphere, because the waves 
are partially reflected and also because of the higher phase 
speed the wave energy spreads out more quickly. In order 
to visualize such waves, we therefore subtract consecutive 
time instances from each other to remove the background 
as much as possible. This essentially amounts to taking 
the time derivative. Thus, in the plots shown below, for any 
quantity X, we show the difference AX = X{t) — X{t — At), 
where At is chosen to be 30 seconds. Figure 1 shows the time 
evolution of the total magnetic field changes (AB), in nT, in 
the noon-midnight meridian plane. The left column shows 
the IOS-1 case, the middle column shows the IOS-2 case, 
and the right column shows the FPS case. Each row shows 
different time, and the time difference between rows is three 
minutes. The red color in the color bar shows an increasing 
magnetic field B, and the blue color indicates where B de¬ 
creases. The color bar range is ±2 nT. In all columns, the 
first plots show the instant right after the FFSs crosses the 
bow shock. 


Table 1. Upstream and downstream plasma parameters for 
the inclined oblique case IOS-1 with shock speed of 380 km/s 
and impact angle 6^^ = 30°. Upstream Mach number and 
obliquity are also shown. 


IOS-1, Va 

= 380 km/s, M 

= 3.7 



0b„ - 

= 51°, 

9x„ = 30' 

D 



Bx 

B, 

Vx 

Vz 

P 

n 

upstream -1.83 

-6.83 

-400.00 

0.00 

20.0 

5.0 

downstream -0.52 

-9.09 

-434.15 - 

■17.65 

67.45 

7.5 


® B in nT, v in km/s, P in pPa, and n in particles/cm®. 


Table 2. Same plasma parameters for the inclined oblique 
case IOS-2 with shock speed of 650 km/s and the same impact 
angle. 


IOS-2, 

Va = 650 km/s, 

M = 7.4 



dB, 

^ = 45°, dxr, = 

30° 



Bx 

Bz Vx 

Vz 

P 

n 

upstream -1.83 

-6.83 -400.00 

0.00 

20.0 

5.0 

downstream -0.52 

-9.09 -461.53 

-28.61 

109.74 

7.5 


Table 3. Same plasma parameters for the frontal perpendic¬ 
ular case FPS with shock speed of 650 km/s and impact angle 
= 0 °. 


FPS, Va 

= 650 

km/s, M 

= 3.7 




= 90° 

o 

II 

s 




Bx 

Bz 

Vx 

Vz 

p 

n 

upstream 0.00 

-7.07 

-400.00 

0.00 

20.0 

5.0 

downstream 0.00 

-10.61 

-483.33 

0.00 

191.37 

7.5 
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Figure 1. 5 consecutive frames representing the resnlt of numerical simulations of AB(nT) plotted in 
the meridian plane with At = 3 minutes, in a time range of 12 minutes. From left, first column shows 
the IOS-1 case, central column represents the IOS-2 case, and the last column shows the FPS case. See 
text for details. 
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Figure 2. Difference of field-aligned currents, AFAC(^tA/m^), for the northern hemisphere ionosphere 
in the same sequence as represented in Figure 1. The center of each plot is the magnetic pole. The right 
side of each plot is dawn while the top is noon (or towards the Sun), and the bottom is the midnight. 
See text for details. 
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Figure 3. DAPEF(mW/m^) plotted in a time range of 20 minutes with time interval of 5 minutes. The 
first column represents the IOS-1 case, second the IOS-2 case, and third the FPS case. Note that only 
the FPS case shows the occurrence of a substorm onset. See text for more details. 
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In Figure 1 the shock fronts appear much broader 
than they really are. First, by taking differences 30 sec¬ 
onds apart, the shock propagates 2-3 Re across the 
grid during that time, thus the shock will appear at 
least that wide. Second, because of numerical diffusion, 
the shocks have a foot or ramp on either side of the 
shock front. Normally, this would be hardly visible in 
a color-coded plot. However, because we take the dif¬ 
ferences and clamp the color bar at low values, these 
numerical artifacts become emphasized and make the 
shock appear much wider than it really is. In the sim¬ 
ulation, the shocks are resolved to within 2-3 cells, i.e., 
to less than 0.5 Re within 10 Re from the x-axis. 

As soon as the shock impacts the magnetopause, it 
launches Alfven waves and magnetosonic waves into the 
magnetosphere. Because the wave speeds are much 
higher in the magnetosphere, these waves race ahead 
of the IP shock in the magnetosheath. The second row 
shows for both cases the time just after the IP shock 
has impacted the magnetopause. In the lOS-I and 
IOS-2 cases, the first contact between the shock and 
the magnetopause occurs just past the northern cusp. 
The induced wave propagates through the northern lobe 
and reaches the plasma sheet in the nightside from the 
north. At this time, the IP shock has not yet impacted 
the southern hemisphere, and consequently there is no 
corresponding wave in the southern hemisphere yet. 

The FPS case (second row, right panel) is distinc¬ 
tively different. The impacts on the northern and south¬ 
ern lobes occur simultaneously, and symmetric waves 
are launched from both the northern and southern mag¬ 
netopause into the magnetosphere. These waves con¬ 
verge on the tail plasma sheet and cause a more signif¬ 
icant compression of the plasma sheet. 

As time progresses, these waves propagate further 
tail-ward. In the asymmetric cases, the waves from 
the southern hemisphere reach the near-Earth plasma 
sheet approximately 3 minutes after the waves from the 
northern hemisphere. As a result of such asymmetric 
impact, there is much less compression of the plasma 
sheet. Instead, the entire plasma sheet is bent south¬ 
ward. This is best visible at the later times (fourth 
and fifth row), where the deflection at x=-30 Re is as 
much as 3 Re in the lOS-I case and 4 Re in the IOS-2 
case. By contrast, in the symmetric case there is no 
such deflection, but instead a transient compression of 
the plasma sheet. 

After the IP shock has passed (bottom row), the mag¬ 
netosphere state seems to be very similar for the three 
cases. However, the particular display only shows dif¬ 
ferences to highlight transients and thus provides little 
information abut the state. 

In order to examine the effects of the IP shocks on 
geomagnetic activity, we examine relevant ionospheric 
quantities. Figure 2 shows the time differences of the 
field-aligned current density (FAC, in /xA/m^) in the 
northern hemisphere polar cap, displayed in the same 


way as the magnetic field evolution of Figure I. The 
times in Figure 2 correspond to the times shown in Fig¬ 
ure I, i.e., a time range of 12 minutes, plotted in three 
minute increments. The range in the color bar is ±0.8 
^A/m^, and regions in red indicate a positive change in 
FAC, while regions in blue indicate a negative change 
in FACs. The dashed circles represent the magnetic 
latitude Am in 10 degree increments from 55° to the 
pole. Left, middle, and right columns represent the re¬ 
sults for the IOS-1, IOS-2, and FPS cases, respectively. 
In all cases, in the first plot, the ionosphere is steady 
because the FFSs have not yet impacted the magne¬ 
topause. Also, in the three situations, the first FAC 
changes are seen about three minutes after the shock 
impact, as can be seen at t = 20:00 minutes for the in¬ 
clined cases and t = 21:30 minutes for the FPS case. In 
the two lOS cases, the FAC changes are mostly in the 
vicinity of the cusp, whereas the FFS causes a broader 
signature that encompasses the entire dayside auroral 
region. The changes in the lOS-I and IOS-2 cases are 
very similar geometrically, with the only difference that 
the IOS-2 response is slightly stronger than the lOS-I 
response. 

As time evolves, with the exception of the first min¬ 
utes after the shock impact, the FAC signature dimin¬ 
ishes in all cases, as can be seen in the middle row. 
Subsequently, activity in the nightside develops for the 
three cases. However, in the FPS case, the ionosphere 
response is much stronger. At t = 24:30 minutes, the 
enhancement of FACs is most evident on the nightside 
of the ionosphere between Am = 70° and Am = 65°, and 
between 0300 MLT and 0600 MLT. At t = 27:30 min¬ 
utes and t = 29:30 minutes, FACs are enhanced close to 
midnight local time, which is a typical substorm signa¬ 
ture [Akasofu, 1964; McPherron, 1991]. The last row of 
Figure 2 shows the strongest nightside FAC variations 
occurring for the FPS case at t = 29:30 minutes, cover¬ 
ing almost all the nightside ionosphere for Am between 
70° and 65°. We attribute this activity due to substorm 
activity that was triggered by the converging waves in 
the plasma sheet. 

Figure 3 shows the diffusive auroral e“ precipitation 
energy flux (DAPEF) in the time steps of 5 minutes 
from t = 25:00 minutes. In the IOS-1 and IOS-2 cases 
the DAPEF remained almost unaffected by the shock 
impact. In all cases, auroral precipitation is enhanced 
in the night side ionosphere approximately 12 minutes 
after shock impact. In the lOS cases, the shocks hit the 
magnetosphere behind the cusp, leading to a mild au¬ 
roral precipitation at t = 35:00. However, at the same 
instant, the FPS enhances more auroral activity in the 
ionosphere nightside because the FPS impacts the mag¬ 
netosphere behind the cusp, and thus the cusp is neither 
displaced nor compressed. In the FPS case, the shock 
hits the magnetosphere first at the nose, leading to a 
compression of the magnetosphere and the cusp, and to 
a pole-ward displacement of the cusp. 



Integrated FAC(MA) Integrated FAC(MA) 
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Later precipitation changes all occur in the night- 
side. As already shown by the FACs, these changes 
are mostly related to substorm activity. In the FPS 
case, auroral substorm onset is formed in the ionosphere 
nightside at 70° < Xm < 65° between 0300 MLT and 
0600 MLT. Such auroral activity is not found in any 
lOS case. 

In order to perform a quantitative comparison we 
integrate the ionosphere quantities over the northern 
hemisphere polar cap. To separate the directly driven 
response, which occurs mostly in the dayside, from the 
induced substorm response, which affects mostly the 
nightside, we integrate the FACs separately for the day- 
side and the nightside. 

Figure 4 shows the time series of the integrated FACs, 
in MA, over the northern hemisphere ionosphere on 
both dayside (a) and nightside (b), for the three cases. 
In both panels, the blue lines indicate the IOS-1 case, 
the green lines indicate the IOS-2 case, and the red 
lines indicate the FPS case. The first vertical dashed 
line indicates the instance at which the IP shock hits 
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Figure 4. Total field-aligned currents (FACs, in MA) in¬ 
tegrated in the northern hemisphere ionosphere. The top 
panel (a) shows the time evolution of the integrated FACs 
in the dayside ionosphere. The bottom panel (b) repre¬ 
sents the integrated FACs in the nightside ionosphere. In 
both plots, the first dashed vertical line at t=16:45 min 
indicates the instant of impact of the IOS-1 and IOS-2. 
The second dashed vertical line, at t=18:28 min, indi¬ 
cates the instant of impact of the FPS. 


the bow shock in the IOS-1 and IOS-2 cases, and the 
second vertical dashed line indicates that instance for 
the FPS case. The plotted quantity is the magnitude of 
the FAC, which is primarily the Region I current. As 
expected, the IP impact enhances the FACs on the day- 
side. Before the shock impact, the FACs on the dayside 
have nearly the same magnitude, 0.29 MA, and were in 
a quasi-steady state. In all cases, the dayside FAC mag¬ 
nitude begins to rise approximately 6 minutes after the 
shock impact. The FACs increase nearly linearly over a 
period of approximately 4 minutes, which corresponds 
roughly to the time it takes for the shock to pass over 
the dayside magnetosphere. This initial rise is nearly 
identical for all cases. Later, in the two lOS cases, the 
FAC remains by and large steady at around 0.5 MA 
for the lOS case, but continues to rise for another 20 
minutes in the FPS case. 

Figure 2 shows that the initial rise is in all cases due 
to enhanced FACs in the vicinity of the cusp. How¬ 
ever, in the IOS-1 and IOS-2 cases the enhancement is 
mostly in the nightside of the cusp, whereas in the FPS 
case FACs both on the dayside and the nightside of the 
cusp are enhanced. It seems that the more thorough 
compression of the dayside magnetosphere leads to the 
stronger and more long lasting FAC enhancement in 
the FPS case, although that is not directly apparent in 
Figure 2. 

The nightside current enhancements (Figure 4 (b)) 
begin about 2 minutes after the dayside enhancements, 
in the three cases. For all shocks, the enhancements are 
qualitatively different from those on the dayside. In all 
cases, strong ULF waves are excited, with a period of 
4-5 minutes. In the two lOS cases, the average FAC 
is enhanced from ~2 MA to ~3 MA due to the shock 
impact, with superimposed waves with an amplitude of 
the order of ^0.2-0.3 MA. In the IOS-2 case, 17-18 min¬ 
utes after the shock impact, the ULF wave amplitude 
rises to ~0.75 MA. In the FPS case, the increase of the 
average FAC is similar, but the wave amplitude is much 
larger, i.e., of the order of ^^1-2 MA. These large oscil¬ 
lations persist for at least 60 minutes. While much of 
the response is similar to what Guo et al. [2005] found 
in their simulations, the oscillations are new. It is at 
present not clear what these oscillations correspond to 
in the magnetosphere, but they are likely cavity modes 
[Samson et al., 1992]. 

Figure 5 shows the cross polar cap potential (CPCP) 
and the integrated precipitation energy DAPEF in pan¬ 
els (a) and (b), respectively. DAPEF is calculated as the 
thermal energy flux of plasma sheet electrons assuming 
perfect pitch angle scattering, and a full loss cone [see 
Raeder et al, 1998; Raeder, 2003, for details].The blue 
line represents the IOS-1 case, the green line represents 
the IOS-2 case, and the red line represents the FPS case. 
Vertical, dashed lines, as described above, indicate the 
shock impact times for all cases. Before the shock im¬ 
pact, the system oscillated noticeably in an amplitude 
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less than 5 kV in both cases. After the impact, the 
IOS-1 and IOS-2 induced a very similar potential jump 
from roughly 35 kV to a peak of approximately 55-60 
kV. The CPCP oscillated with a period of nearly 5 min¬ 
utes until it reached the near steady state value of 25 
kV after 42 minutes. 

The FPS case is more dynamic. The system oscil¬ 
lated around 30 kV before the shock impact. Right two 
minutes after the FPS hit the bow shock, the poten¬ 
tial difference dropped to 23 kV. This effect was also 
seen by Guo et al. [2005] and interpreted as the redis¬ 
tribution of the FACs after the shock penetrates the 
magnetosphere. The potential drop then reached two 
smaller peaks, 50 kV at t=22:00 minutes, and 62 kV 
at t=24:00 minutes until it reached the maximum of 
82 kV at t=27:00 minutes. The potential still reached 
two peaks of 53 kV at t=30:00 minutes and t=33:00 
minutes. Then, the potential difference decreased in 
a period of nearly 7 minutes. After close to t=33:00 
minutes, both systems evolved to nearly the same final 
quasi-steady state. This similarity has been shown by 
Guo et al. [2005] but not with this oscillatory behavior. 

The bottom panel (b) of Figure 5 shows the time se¬ 
ries for the integrated DAPEF, in units of GW. Again, 
the blue line represents the IOS-1 case, the green line 
represents the IOS-2 case, and the red line represents 
the FPS case. Before the shock impact the ionosphere 
is in a quiet state with low auroral activity. After the 
shock hit the bow shock, between t=24:00 minutes and 
t=33:00 minutes, the three systems evolved quite simi¬ 
larly, while the auroral energy flux suffered a small drop 
in the FPS case at t=34:00 minutes. In the lOS cases, 
the DAPEF attains a maximum value of barely 100-120 
GW near t=44:00 minutes, and then evolves to a final 
state of ~20-50 GW. On the other hand, the FPS en¬ 
hanced the DAPEF much more, and briefly reaches a 
maximum of ~350 GW at t=45:00 minutes. 

By comparing the two shocks with the same Mach 
number, namely IOS-1 and FPS, we observe that both 
shocks lead to very different geomagnetic responses. 
Also, the IOS-2 geomagnetic response is smaller in com¬ 
parison to the FPS geomagnetic response, even though 
the inclined shock was twice stronger than the head-on 
shock. We attribute these results to the different shock 
normal orientations. Gomparison with Figure 3 shows 
that this enhanced precipitation flux must be from the 
nightside due to substorm activity. Such precipitation 
enhancements have been reported earlier by Zhou and 
Tsurutani [2001], Tsurutani and Zhou [2003], and Yue 
et al. [2010], who find that the auroral precipitation and 
field-aligned currents in the nightside can be intensified 
after a FFS impact, because it triggers the release of 
stored magnetospheric/magnetotail energy in the form 
of particularly large substorms, or even supersubstorms 
[See Tsurutani and Lakhina, 2014, for more references 
therein]. Such substorm triggering might be a result 
of the decreasing in the nightside Bz: as observed at 
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Figure 5. (a) Cross polar cap potential (CPCP, in kV) 
in the northern hemisphere ionosphere, (b) Integrated 
auroral electron precipitation energy flux (DAPEF, in 
GW) in the same hemisphere. The dashed vertical lines 
are the same as in Figure 4. 

geosynchronous orbit by Sun et al. [2011, 2012]. In 
that case, the decrease in the nightside Bz is suggested 
as a result of Earth-ward transportation of magnetic 
flux by temporarily enhanced plasma flows at the night¬ 
side geosynchronous orbit after the impingement of an 
IP shock. In our case, the symmetric plasma trans¬ 
port leads to a more intense geomagnetic activity. Al¬ 
though our simulations do not represent a storm and 
although the simulation results may be qualitatively in 
error, they match qualitatively the observed shock im¬ 
pact behavior. 

4. Summary and Conclusions 

It has been known for a long time that IP shocks can 
have a profound impact on the magnetosphere; how¬ 
ever, it is much less known which factors determine the 
geoeffectiveness of IP shocks. There have been a few 
studies in the past addressing the effect of shock geom¬ 
etry [Takeuchi et al., 2002; Jurac et al., 2002; Guo et al., 
2005; Wang et al., 2005, 2006; Grib and Pushkar, 2006; 
Samsonov, 2011], but none have considered the partic¬ 
ular geometry that we investigated here. Specifically, 
here we use global simulations to contrast two differ¬ 
ent scenarios. In one scenario, the IP shock normal lies 
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along the Sun-Earth line, such that there is a frontal 
impact on the magnetosphere. In the other scenario, 
the IP shock is inclined with respect to the Sun-Earth 
line. In either case, the shock normal lies in the GSE 
x-z plane, and the IMF is southward, such that there is 
no y-dependence of any solar wind parameter. The two 
scenarios lead to very different responses of the magne¬ 
tosphere: 

1. In the frontal case, the shock launches waves sym¬ 
metrically into the magnetosphere, which converge on 
the tail plasma sheet and compress it. By contrast, in 
the inclined cases, the waves reach the plasma sheet 
at different times, causing much less compression, but 
a north-south displacement of the plasma sheet. This 
result holds even for shocks with larger Mach numbers. 

2. In the frontal case, the compression triggers a sub¬ 
storm, whereas in the inclined cases there is no excess 
geomagnetic activity beyond what would be expected 
for southward IMF. 

3. In all cases the shock impact enhances FACs in 
the dayside with a similar quantitative response. How¬ 
ever, in the inclined cases the FAC enhancement occurs 
mainly behind the cusp, whereas in the frontal case the 
cusp is displaced while the FACs increase over a wider 
MET range. In the frontal case, the dayside FAC re¬ 
sponse also persists longer and is more intense, i.e., a 
~200% enhancement versus a ~I00% enhancement in 
the inclined case. 

4. The nightside FAC response is qualitatively sim¬ 
ilar in the three cases and shows the development of 
ULF waves. However, in the frontal case the ULF wave 
amplitude is much stronger. The detailed excitation 
mechanism remains to be investigated. 

5. The response of the cross polar cap potential is rel¬ 
atively benign and limited to the first 15 minutes after 
the impact, in all cases. The three cases relax to the 
pre-impact state in less than 20 minutes. 

6. In all cases diffuse auroral electron precipitation 
increases in similar fashion in direct response to the 
shock impact. In the frontal case, this is followed 
by a delayed response ~I0 minutes later, which peaks 
~20 minutes after the impact, and which comes from 
the nightside. The latter is interpreted as a consequence 
of the substorm triggered by the shock impact. 

Our results show that the shock impact angle has a 
major effect on the geoeffectiveness of the shock, even 
more than the Mach number or some other measure of 
the shock strength. Although we only covered a rela¬ 
tively small parameter space in terms of impact angles, 
shock strength, and IMF orientation, the qualitative 
and quantitative differences we found are signihcant. 
With respect to substorm triggering [Kokuhun et al., 
1977; Lyons, 1995, 1996; Lui et al., 1990], the differ¬ 
ences we found in the inclined and frontal cases are of 
particular importance. Apparently, the same type of 
IP shock can either trigger a substorm or not, depend¬ 
ing on the shock normal direction, which has not been 


considered in previous studies. We also find large am¬ 
plitude waves in the FACs that are apparently caused 
by the shock impacts. Civen their period, these waves 
are likely cavity modes [Samson et al., 1992]. We find 
that the modes have significantly larger amplitude for 
the frontal case. This is likely due to the fact that 
the waves that converge on the tail and compress the 
plasma sheet from there launch a wave back towards 
the nightside magnetosphere, which in turn excites the 
cavity mode. In the inclined cases, this earthward wave 
is likely much weaker, and thus excites a weaker cavity 
mode. We will study the wave excitation in more de¬ 
tail in a forthcoming paper. Other future work is also 
clearly laid out, namely finding the correlation between 
geomagnetic activity and IP shock normal orientation 
in data, and a better parameter space coverage with 
simulations. 
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